Self-duality in quantum impurity problems 
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We establish the existence of an exact non-perturbative self-duality in a variety of quantum impurity 
problems, including the Luttinger liquid or quantum wire with impurity. The former is realized in 
the fractional quantum Hall effect, where the duality interchanges electrons with Laughlin quasipar- 
ticles. We discuss the mathematical structure underlying this property, which bears an intriguing 
resemblance with the work of Seiberg and Witten on supersymmetric non-abelian gauge theory. 
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Duality is a rather ancient concept. It was probably 
observed first in electromagnetism, where the theory is 
invariant under exchange of electricity with magnetism. 
For the semi-classical theory to remain local, the electric 
charge e and the magnetic charge m obey the Dirac quan- 
tization condition em = 27r, so a large e (which plays the 
role of the coupling) can be traded for a small m. In 
general, duality maps a theory with strong coupling to 
one with weak coupling, and is thus a powerful tool for 
exploring strongly-coupled regimes. 

A theory is self-dual when there is an exact map from 
a theory at strong coupling to the same theory at weak 
coupling. The best-known example occurs in the two- 
dimensional Ising model, where the map is known as 
Kramers- Wannier duality, and provides a quick way to 
obtain the critical temperature. 

In this paper, we present a exact form of self-duality in 
1 + 1 dimensional quantum impurity problems. This self- 
duality is much more powerful than in the Ising model, 
because there is no intermediate-coupling critical point 
interfering with the analytic continuation from weak to 
strong coupling. This duality had sometimes been used 
in the past as "approximate" , with no clear understand- 
ing of its exact status. We show here that it results from 
integrability, a powerful symmetry which restricts drasti- 
cally the form of the Lagrangian in the low energy limit. 
We make the duality very explicit by obtaining simple 
integral expressions for transport properties. These hint 
at a deep mathematical structure, somewhat similar to 
the one arising in supersymmetric gauge theories with 
Seiberg- Witten duality 0. 

We study a model of N flavors of interacting gapless 
fermions moving in one dimension. We denote the left- 
moving and right-moving modes by ipLi(x) and ipm{x) 
respectively with i = 1 . . . N. We consider an impurity 
located at x = which gives a backscattering interaction 
Lb = Xip Li (0)ipRi(0)+cc. These models can be bosonized 
in terms of N free bosons, so that the only interaction 
between left- and right-moving modes is at x = 0. For 
N=l, this is the Luttinger liquid with impurity, which 
in its bosonic formulation is usually called the boundary 



sine-Gordon model. For N—2, this is a quantum wire; 
the two species correspond to the two spins up and down 
of the electrons. We consider only the case where the 
N species are identical, and we allow for charge inter- 
actions, parametrized by a coupling constant g (g — N 
at the free point). Since the excitations are gapless in 
the bulk, the only scale in the problem is the coupling 
A from the impurity interaction. There are two critical 
points: the "no-backscattering" point where A = 0, and 
the "strong-coupling" fixed point where A — > oo and all 
the excitations are backscattered. 

Let us first discuss the simplest case N=l, the Lut- 
tinger liquid with impurity. As is well known, the inter- 
acting fermions of the Luttinger model can be bosonized 
in terms of a single free boson. Moreover, one can "fold" 
space in half around the impurity at x—0 and restrict 
to the half- line < x < oo Without the impurity, 
the Lagrangian is L = (l/4ng) f^°dx (d^cf)) 2 where g 
parameterizes the interactions of the original Luttinger 
fermions. The two critical points correspond to Neumann 
boundary conditions d x (j)(0,t) = and Dirichlet bound- 
ary conditions </>(0, t) = const. As is familiar from confor- 
mal field theory, the massless boson can be decomposed 
into left- and right- moving pieces <f> = <f>L + 4>R- One 
can then define the "dual" boson <f> — (<f>L — <t>R)/g and 
rewrite the Lagrangian as L = (g/An) J °° dx (<9 M 0) 2 . 
Dirichlet boundary conditions on <f> correspond to Neu- 
mann boundary conditions on (f>, and vice versa. With no 
interactions at the boundary, the theory is therefore in- 
variant under the mapping g — > 1/g, and the interchange 
of Dirichlet and Neumann boundary conditions. 

An interesting set of charge-carrying operators are 
written in terms of this boson and its dual as ^ e . m (x, t) = 
exp[ie(j)(x, t) + im<f>(x,t)]. The parameters e and m 
must be integers to ensure locality. For Euclidean time 
this entire problem can be reformulated as a classi- 
cal two-dimensional Coulomb gas H, where "i+o cre- 
ates an electrical charge e, and ^0,1™ creates a mag- 
netic charge m. In this context, the condition that e 
and m be integers is analogous to the Dirac quanti- 
zation condition. The boundary two-point function is 
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(*"e,m(0, 0)*_ e ,— m(0, t)) ~ t~ 2a!e > m where the boundary 
scaling dimension x e . m = (e^/g + m/^/g) , which has the 
electric-magnetic duality g — > l/g. 

Including the impurity interaction causes operators to 
be added to the Lagrangian. A relevant boundary oper- 
ator has boundary dimension x e _ m < 1. Without loss of 
generality we can consider g < 1, because any g > 1 
model can be mapped to a g < 1 model by sending 
g — > l/g and exchanging the two fixed points. This du- 
ality is not self-duality; it does not say anything about 
the g < 1 models by themselves but only that the g > 1 
models are equivalent. When g < 1, only operators with 
to = are relevant [||. When g > 1/4, only ^fifi and its 
conjugate VP-i.o are relevant. Therefore, we consider the 
backscattering Lagrangian 

L B = A(*!, Q + *_i,o) - 2Acos0(O,*). (1) 

This interaction induces a flow from Neumann bound- 
ary conditions (UV) to Dirichlet (IR) as A increases, 
or, equivalently, as energy decreases. From this expres- 
sion one finds for weak backscattering at zero tempera- 
ture, the current / and voltage V obey the scaling law 
I-gV oc V 29 - 1 §). 

The self-duality appears in the action in the strong 
coupling limit. At very large A, the field 0(0, t) is pinned 
at the minima of the potential, so it is reasonable to study 
the vicinity of the Dirichlet fixed point by considering 
the instantons interpolating between neighboring minima 
|^| . After some elementary computations, one obtains the 
leading (charged) irrelevant operator 

L B « Ad(* ,i + *o,-i) = 2Az,cos0(O,i). (2) 

The subscript D is for dual, because this form of the 
backscattering lagrangian is dual to the form ([!]). Since 
cos 0(0, t) has dimension l/g, Xo oc A -1 / 9 , and at zero 
temperature the current and voltage obey the scaling law 
/ cx V 2 ^ 9 ^ 1 in the strong-backscattering limit Q]. Thus 
the strong and weak-backscattering exponents are related 
by interchanging g — ► l/g. 

If (||) were exact, the theory would be completely self- 
dual under the interchange A — > Xn and g — ► l/g. This 
has been often assumed in the literature, and is of course 
very appealing, especially in the context of the fractional 
quantum Hall effect at v = 1/3. Here, the Luttinger 
fermions at g — 1/3 are equivalent to Laughlin quasipar- 
ticles, and the impurity corresponds to a point contact 
allowing tunneling between edges of a fractional quan- 
tum Hall device ||. In the weak-backscattering limit, 
the tunneling ([!]) is caused by Laughlin quasiparticles. 
In the strong backscattering limit the system is effec- 
tively two different Hall devices, weakly coupled by (||) 
so that only electrons can tunnel. Thus the duality maps 
a description in terms of Laughlin quasiparticles to one 
in terms of electrons. 



However, there is no reason a priori why the duality 
should be exact: the instanton expansion of is not 
exact, and multi-instanton processes can be expected to 
contribute, adding to (||) terms oc cosn0(O,f). In addi- 
tion, neutral operators will also contribute, for instance 
the density - density coupling (50) 2 (which has dimension 
2 and so for g < 1/2 is the leading irrelevant operator). 
The perturbative computations near the strong coupling 
fixed point will involve all of these terms, and thus have 
little to do with its weak-coupling dual (|l|). 

It turns out however that Bethe ansatz computations 
established that there is indeed an exact duality for 
the current at zero temperature: 

I(\ 1 g 1 V)=gV-gl(\ D ,~,gv\ (3) 

Some well-motivated and well-checked conjectures re- 
quire that this also hold at arbitrary T 0. The self- 
duality relation (^) can be derived directly by imposing 
the strong constraints of integrability on the IR action. 
The integrability of the boundary sine-Gordon model 
means that there are an infinite set of conserved quanti- 
ties in involution, even with the perturbation (Q) of the 
UV fixed point. These are generated by the non-local 
operator cos0, and by an infinite set of local operators 
02k+2 which are, roughly, powers of the stress energy 
tensor, with dimensions 2fc + 2, k an integer. By consid- 
ering the Yang-Baxter equation and its various algebraic 
consequences in detail [|7|]lC|] , or more simply by demand- 
ing integrability of the flow both near the UV and near 
the IR fixed point, it can be shown that the complete IR 
lagrangian is given by (^) plus an infinite series of the 

tyP e YlkLo Ck^D +r>9 X 0%k+2 (this sort of statement 
of course can only make sense within a particular regu- 
larization scheme, here dimensional regularization ||10f ). 
The remarkable fact is that only one charged operator 
appears. It is also crucial that the neutral operators all 
commute: as a result, it can be shown, by using the non 
equilibrium Keldysh formalism, that DC transport prop- 
erties are insensitive to anything but (||), and thus are 
self-dual, at any temperature. Therefore (^) must hold 
to all orders in perturbation theory. 

The existence of the duality prompts questions about 
the analytic structure of the whole problem: knowing a 
priori that (|J) holds, could one find the exact formula 
for I without using the Bethe ansatz? We present here 
a formulation of the problem that suggests it should be 
the case, and bears an intriguing resemblance with M 
(for a previous attempt in a different direction, see ||). 
Introducing the dimensionless parameter u oc VX 1 ^ 9 ^ 1 ^ 
(the precise relation is given in &), doing perturbation 
theory in powers of A in (|l|) gives 

J- = 1-J2a n (g)u 2 ^-V (4) 
^ n=l 
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for large u (small A). The analogous expansion in the 
strong-backscattering regime (small u) must be 
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2n(l/ S -l) 



(5) 



because the self-duality relation (H) means that not only 
the powers but also the coefficients in the expansions of 
the current are dual under g — > Considering the di- 
mensionless current T(g,u) = I/gV as an analytic func- 
tion of gives the non-perturbative self-duality 



1{g,u) = 1 -I(l/g,u). 



(6) 



This equation applies for all g, so does much more than 
just relate a model at g > 1 to one with g < 1. 

Given either one of the expansions @) or (||), we can 
now write an expression for the current which not only 
makes proving the self-duality simple, but which displays 
an intriguing connection to Seiberg-Witten theory. The 
Bethe ansatz gives the coefficients a n to be 



a n (g) 



(-1) 



n+l 



T(ng + 1) 



2 r(3/2 + n(g - 1)) 



(7) 



With this a n , the expansions (J|) and (]5|) have finite ra- 
dius of convergence, so the self-duality must hold non- 
perturbatively. Simple expressions for the magnetization 
can be found by using the integral form of the B func- 
tion r(a)r(6)/r(a + b) to rewrite a n (g). Since the sum 
(|J) is absolutely convergent for large enough u, we can 
interchange the order of summation and integration, and 
do the sum. This yields our result 



Z(g,u) = — I dx 
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\J x + x 9 



(8) 



where the contour of integration Co starts at the origin, 
loops around the square-root branch point on the pos- 
itive real axis, and returns to the origin. There must 
be such a branch point because the expression inside 
the square root is positive for large x and negative at 
x = 0. By using the fluctuation-dissipation theorem 
((A/) 2 ) = {{gV) 2 /{l-g))d v l derived in §, one obtains 
an integral expression for the DC shot noise as well. 

It can easily be checked that the series (jj) can be re- 
covered from (^) by changing variables in the integral to 
t = x/u 2 , and expanding the resulting integrand in pow- 
ers of u 2g ~ 2 . However, the integral in (||) converges for all 
real u, not just in the large- w regime where the sum (^) 
does. Therefore it must also yield the expansion (pft in the 
small- it limit. By defining the variable t — xjv?! 9 and 
expanding the resulting integrand in powers of u 2 / g ~ 2 , 
one can indeed verify this fact. Moreover, one can prove 
the self-duality relation (^) without the series expansions, 
because the integral expression (^) is analytic in g and 
so can be utilized for any g: 
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dx- — . 

Co V X + X 9 



dx-^- 

ox 



= 1 -l{g,u) 



Co 



\/ X + X 9 — 



l(g,u) 



(9) 



where to get from the first to second line we change vari- 
ables x — > x g , and Cq changes accordingly. 

Our expressions for the current bear a strong resem- 
blance with the representations obtained in (!]] for quan- 
tities such as mass gap in N=2 supersymmetric gauge 
theory. There the parameter u is an order parameter re- 
lated to the expectation value of the Higgs field. This 
is perhaps not so shocking. The idea of Seiberg-Witten 
theory is to exploit the fact that certain physical quan- 
tities are analytic in some coupling, except at specific 
singularities. If one takes the coupling around a singular- 
ity, physical quantities do not necessarily return to their 
original values (for example if they have logarithmic be- 
havior near the singularity): there is a non trivial Berry's 
phase, or monodromy. If one knows all the singularities 
and monodromies of an otherwise-analytic function, one 
can reconstruct the complete function for any value of 
the coupling: its typical form in the present context is 
given by J dx/y, where y 2 is a rational function of x, and 
the integral loops around some of the square-root branch 
points of the integrand. The coefficients in y 2 depend on 
the couplings of the theory. Singularities occur at val- 
ues of the couplings where two roots of y 2 coincide, so 
that the integral logarithmically diverges. In the gauge 
theories discussed in Q and in many subsequent gener- 
alizations, one can use perturbation theory to find the 
monodromies. Thus perturbative information and ana- 
lyticity leads to non-perturbative results. 

In integrable models, analyticity in couplings has long 
been exploited to do various computations |p.l| . In the 
present context, the analyticity is in the dependence 
of the current on complex u (that is, voltage). Little 
seems to be known in general about this question, but 
the integral expression (||) reveals a very simple struc- 
ture. The current is singular at the values uq where 

0, namely 

4 = (- 

of: 



both x + x 9 — Uq = and I + gx 9 
■ 2 = (-n\a/{l-9){i - 5 ). when u = u , two of the roots 
x + x a — u 2 coalesce. The integral (||) diverges if the 
contour runs in between the two coalescing roots; this 
never happens for positive real u. Thus we know the 
exact location of all of the Lee- Yang zeroes of the cur- 
rent. All of them have the same absolute value \uq\, and 
for rational values of g, there are a finite number (for 
some reasons why this is, see Q). The value |uo| gives 
the radii of convergence: the series (Q) absolutely con- 
verges for |u| > |uq|, while (||) absolutely converges for 
\u\ < \u \. 
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There are also non-trivial monodromies in this prob- 
lem that can be extracted from the integral expressions. 
For example, when g = 1/2, the integral can be done to 
find the known result T = 1/2 — (l/2u) arctan(2u). This 
has singularities at u$ = ±i/2. Looping u around one of 
the singularities moves onto a different sheet of the arc- 
tangent, so I — > X ± ir/(2u) for g = 1/2. Unfortunately, 
for general g the monodromies are more complicated and 
at the present time we do not know how to obtain them 
without already knowing the integral (||). 

The foregoing self-duality extends to the case of N 
channels. The case N=2, a quantum wire with impu- 
rity, was studied perturbatively in Q , and the exact cur- 
rent was found using the Bethe ansatz in jl3| . Following 
|fl^3|, it is straightforward to find the large- m series ex- 
pansion of the current for general N. Defining as before 
2n = I/gV, g the conductance without impurity (g = N 
for N channels of free electrons), one finds 



z"w = £(-i) 



71=0 



v^Fr(n(l + Nh) + l)T(nh + l) 
~2~ n\r{nNh+l)T(nh + 3/2) ' 



,2nh 



where h — - — -Jr. This can be written as the integral 

g N ° 



^n(9, u) — ~r \ dx 



{l + x h ) 



h\N-l 



y/x(l+X h ) N -- 



(10) 



This normalized current obeys the self-duality relation 
(^) as with the N—l case, here with the dual coupling 
N 2 /g. It would be quite interesting to see whether this 
result generalizes to the models discussed in |TJ] . 

Duality considerations extend beyond transport prop- 
erties; for instance, it is known from the exact solution 
that there was a sort of duality between the UV and IR 
in the Kondo model |l5[| . This follows the same pat- 
tern as before. Consider the anisotropic Kondo problem, 
which is defined by the boundary Lagrangian general- 
izing d), L B = A(*i )0 »S7 + 9-i t0 S+), where S ± are 



su(2) q spins (q 



— ping 



, g parameterizing the anisotropy) 



in the representation of spin j. As with the tunneling 
problem, the integrability requires the IR Lagrangian to 
be very simple: it contains only the conserved quanti- 
ties 2 k+2 and a single non-trivial term like (|j): L B w 
Xd (4*0. i-Sjli/2 + ^o.-iS^Li/a)) wnere now the spins are 
in the representation of spin j — 1/2. As a result, the 
free energy at magnetic field H satisfies a duality rela- 
tion similar to (@) 



1 1 H 

f(j,x,g,H) = f [j - -,x D , -, — 

2 g g 



The left-hand side is the expansion in the UV, whereas 
the right-hand side is the IR expansion; the extra terms 
on the right-hand side are only in odd-integer powers of 
H/X 1 ^- 1 , and arise because the free energy does couple 
to the charge-neutral operators in the IR action. The 



free energies for all S are also given by integrals like jq) , 
which surround branch points in the complex plane [ fl2[ . 

There are some interesting distinctions between these 
quantum-impurity problems and the Seiberg-Witten the- 
ories. In the gauge theories, couplings never seem to ap- 
pear in the exponents of x. Here g rational leads to a 
polynomial for y 2 as in the gauge theories, but g irra- 
tional is a slightly new situation. However, this is quite 
interesting because none of the g — > 1/g dualities known 
in gauge theory have a proof like that in (|]) . Moreover, 
the Luttinger liquid with impurity is not supersymmetric; 
as argued above the analyticity seems to follow from the 
integrability, not from supersymmetry as in the gauge 
theories. This opens up a new realm of possibilities in 
gauge theory, a few of which are discussed in Q . 

In conclusion, it is likely that duality is a key feature 
of integrable field theories in 1+1 dimensions. For exam- 
ple, a situation similar to what we have discussed should 
occur in the massless flows between minimal models of 
conformal field theory, where we expect a duality be- 
tween the 0i3 perturbations in the ultraviolet and </>3i 
perturbations in the infrared. 
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